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Abstract — This paper studies n-person simultaneous-move 
games with linear best response function, where individuals 
interact within a given network structure. This class of games 
have been used to model various settings, such as, public goods, 
belief formation, peer effects, and oligopoly. The purpose of 
this paper is to study the effect of the network structure on 
Nash equilibrium outcomes of this class of games. Bramoulle 
et al. derived conditions for uniqueness and stability of a Nash 
equilibrium in terms of the smallest eigenvalue of the adjacency 
matrix representing the network of interactions. Motivated by 
this result, we study how local structural properties of the 
network of interactions affect this eigenvalue, influencing game 
equilibria. In particular, we use algebraic graph theory and 
convex optimization to derive new bounds on the smallest 
eigenvalue in terms of the distribution of degrees, cycles, and 
other relevant substructures. We illustrate our results with 
numerical simulations involving online social networks. 

I. Introduction 

In most social and economic settings, individuals do not 
interact uniformly with the rest of a society. Instead, they 
influence each other according to a structured network of 
interactions. This network can represent friendships in a 
social network, transactions among firms in a market, or 
communication links in a process of belief formation. In this 
context, an interesting question is to study how the network 
structure affects the outcome of the interactions of the agents. 
With this purpose, one can model strategic interactions in 
a networked society as a multi-player simultaneous-move 
game. In particular, we focus our attention on the broad class 
of games with linear best response function [1]. This class 
of games have been used to model various settings such as 
belief formation [2], peer effects [3], and public goods [4]. 

In order to analyze the influence of the network struc- 
ture on the game outcome, we use two recent results by 
Bramoulle et al. [5] and Ballester et al. [6] relating the Nash 
equilibria of the game with the largest and smallest eigen- 
values of the adjacency matrix of interactions. For example, 
in [5], the authors show that uniqueness and stability of a 
Nash equilibrium in games with linear best responses can 
be determined by the smallest eigenvalue of the network. In 
[5], it was also illustrated how the smallest eigenvalue of the 
adjacency matrix determines the capacity of the network to 
absorb perturbations on the actions of the agents. 

Motivated by the results of Bramoulle et al. [5], we study 
how local structural properties of the network affect the 
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smallest eigenvalue of the adjacency matrix of interactions, 
affecting game equilibria. Therefore, our results build a 
bridge between structural properties of a network of interac- 
tions and the outcome of games with linear best responses. 
In particular, we use algebraic graph theory and convex 
optimization to derive bounds on the smallest eigenvalue of 
the adjacency matrix in terms of the distribution of degrees, 
cycles, and other important substructures. As we illustrate 
with numerical simulations in online social networks, these 
bounds can be used to estimate the effect of structural 
perturbations on the smallest eigenvalue. 

The paper is organized as follows. In the next subsec- 
tion, we review graph-theoretical terminology needed in 
our derivations. In Section |llj we review the relationship 
between the equilibria of games with linear best responses 
in a network and the smallest eigenvalue of the adjacency 



matrix of interactions. In Section III we use algebraic graph 



theory to derive closed-form expressions for the so-called 
spectral moments of a network in terms of local structural 



features. In Section IV we use convex optimization to derive 
optimal bounds on the smallest (and largest) eigenvalue of 
the interaction network from these moments. Our bounds 
help us to understand how structural properties of a network 
impact the stability properties of the Nash equilibria in the 
game. We illustrate our results with numerical simulations 
in real online social networks in Section [V] 

A. Notation 

Let Q = (V,£) denote an undirected graph with n 
nodes, e edges, and no self-loop^] We denote by V (G) — 
{vi,..., v n } the set of nodes and by £{G)QV (G) x V (</) 
the set of undirected edges of G- If {vi,Vj} E £ (G) we 
call nodes Vi and Vj adjacent (or neighbors), which we 
denote by vi ~ Vj and define the set of neighbors of Vi 
as Mi = {w G V (G) ■ {vi,w} e £ (G)}- The number of 
neighbors of Vi is called the degree of the node, denoted 
by di. We define a walk of length k from vq to Vk to 
be an ordered sequence of nodes {vq,v\, ...,Ufe) such that 
Vi ~ Vi+\ for i = 0, 1, k — 1. If vq = Vk, then the walk 
is closed. A closed walk with no repeated nodes (with the 
exception of the first and last nodes) is called a cycle. For 
example, triangles, quadrangles and pentagons are cycles of 
length three, four, and five, respectively. 

Graphs can be algebraically represented via matrices. The 
adjacency matrix of an undirected graph G, denoted by 
Ag = [a-ij], is an n x n symmetric matrix defined entry- 
wise as aij = 1 if nodes vi and Vj are adjacent, and 

'An undirected graph with no self-loops is also called a simple graph. 



dy = otherwis^] The eigenvalues of Ag, denoted by Xi > 
A2 > . . . > A n , play a key role in our paper. The spectral 
radius of Ag, denoted by p(Ag), is the maximum among 
the magnitudes of its eigenvalues. Since Ag is a symmetric 
matrix with nonnegative entries, all its eigenvalues are real 
and the spectral radius is equal to the largest eigenvalue, Ai. 
We define the fc-th spectral moment of the adjacency matrix 
Ag as: 



m k (Ag) 



1 



(i) 



As we shall show in Section [TIT] there is a direct connection 
between the spectral moments and the presence of certain 
substructures in the graph, such as cycles of length k. 

II. Strategic Interactions in Networks 

In this section we present the game-theoretical model of 
strategic interactions considered in this paper and present 
interesting connections between the Nash equilibria and the 
eigenvalues of the adjacency matrix of the network. 

A. The Model 

We represent the network of influences using a simple 
graph Q. Let M — {1, • • • , n} denote the set of n players 
located at each node of the graph Q. We denote by Xi € 
[0, 00) the action chosen by agent i, and by x the vector that 
represents the joint actions for all agents. We denote by x_i 
the vector of actions of all players excluding player i. As 
mentioned before, players interact according to a network of 
influences that we describe using its adjacency matrix Ag. 
The interactions are assumed to be symmetric, = a 3 », and 
we do not allow self-loops, an = 0. The payoff function for 
agent i is given by: 

Ui(xi,x-i;5, Ag) 

where S € M is a parameter that can be tuned to change the 
influence of neighboring nodes on each player's action. 

B. Games with Linear Best Response Functions 

We study a class of games whose best response functions 
take a linear form. One well known example of this class of 
games is the differentiated-product Cournot oligopoly with 
linear inverse demand and constant marginal cost with payoff 
function defined as [5]: 

t 7 \ 

Ui(xi,X-i; 5, Ag) =Xi I a - b(xi + 26^ a tJ Xj) - dx t , 

(2) 

where d is the constant marginal cost, and xi represents 
the amount produced by agent i in the oligopoly. Here, the 
inverse demand for agent i is given by Pj(xj, x_,; S, Ag) = 
a — b(xi + 2SY^j—i a ij x j)- O ne can prove that the best 
response function for this type of games yield the form [5]: 



/i(x, 5,Ag) = max(0, x t 



"For simple graphs, an = for all i. 



3=1 



(3) 



where x.i is the action that agent i would take in isolation, 
i.e., Xi € argmaxa;. Ui(xi, x_^; 5, Ag) with = 0, for all 
j. Without loss of generality, one can normalize ±i = 1 for 
all i, so that /-j(x) G [0, 1]. Then, a Nash equilibrium for this 
game is a vector x e [0, 1]" that satisfies Xi — /i(x, S, Ag), 
for all agents i G JV, simultaneously. In what follows, we 
briefly describe a strategy to compute the complete set of 
Nash equilibria for 5 E [0, 1]. 

C. Complete Set of Nash Equilibria 

Using the best response function in (J3J, one can determine 
the entire set of equilibria by simultaneously solving for the 
best response of each player. In [5], an algorithm that finds 
the full set of Nash equilibria in exponential time is proposed. 
For a vector x, let S denote the set of active agents, i.e., 
S = {i : Xi > 0}. Let X5 denote the vector of actions of the 
agents in S. The set of active players induce a subgraph Qs C 
Q, with node-set S C V (Q) and a set of edges £ $ C £ (Q) 
connecting active agents. We denote by Qjsf\s,s me subgraph 
of Q whose edges connect active agents in S to inactive 
agents in J\f\S. The adjacency matrices of Qs and G^\s,s 
are denoted by A s € Rl s l x l s l and A^-^s G R"-I S M^, 
respectively. Then, one can show the following [5]: 

Proposition 1: A profile x with active agents S is a Nash 
equilibrium if and only if: 

(J| S | +5A s )x s = 1| S |, and 

$ A JV\S,S x S > l n -|S|) 

where I p is the p x p identity matrix and l q is the q- 
dimensional vector of ones. 

Thus, in order to determine the complete set of all Nash 
equilibria, one can check the conditions in Proposition [TJ for 
each one of the 2™ possibilities of S. For each possible 
S, these conditions can be checked by computing X5 = 
(i|s| + M5) _1 l|5|Qand checking whether SAj^\ S ,S X S > 
1„_|5|. If the last inequality holds, then xg is an equilibrium 
outcome. Note that using this approach to compute the set of 
equilibria runs in exponential time. However, we can analyze 
some properties of the Nash equilibria, such as uniqueness 
and stability, by looking into the eigenvalues of the adjacency 
matrix. 

D. The Shape of Nash Equilibria 

In order to relate the equilibrium outcomes of the game to 
the network structure, the authors in [5] defined the following 
potential function: 



71 / 1 \ 1 71 71 

(x;S,Ag) = ^ [Xi- -x\ ) - J^E fl 



i=l j=l 



Then, they proved, using Kuhn-Tucker conditions, that the 
set of Nash equilibria coincides with the critical points of 

'The matrix (l\s\ + ' s invertible for almost any 8 £ R, excepting 

the measure zero set {— i = 1,...,\S\}, where /Ltigi are the 

eigenvalues of As- 



the following optimization problem: 

(P) max x <p(x;6,Ag) 

s.t. Xi > 0, for all i, 

for a given network structure Ag and a parameter 6. 

E. Eigenvalues and Nash Equilibria 

We can find several results in the literature providing 
sufficient conditions for the existence of a unique Nash 
equilibrium in games with linear best response functions in 
terms of the eigenvalues of the network of influences. We 
enumerate below some sufficient conditions that are related 
with our work: 

Proposition 2: Consider the class of games with linear 
best response functions described in Section |II-B| For these 
games, we have the following sufficient conditions for the 
existence of a unique Nash equilibrium: 

(0 S < -l/X n (A g ), (Bramoulle et al., [5]). 

(if) 8 <l/p (Ag), (Ballester et al., [6]). 

We can compare the set of spectral conditions in Propo- 
sition [2] using the following inequalities [5]: 

Lemma 3: For any simple graph Q, we have that 
— l/X n (Ag) > 1/p (Ag), where this inequality is strict when 
no component of Q is bipartite. 

Remark 4: Hence, Condition (f) in Proposition [2]provides 
the best sufficient condition for the uniqueness of Nash 
equilibrium in these games. Furthermore, one can also prove 
that under Condition (i) or (ii), all players are active at the 
equilibrium point, i.e., S — N [5]. 

(Ag) in 



In Section 



IV 



we shall derive upper bounds on A r , 
terms of structural properties of the network. These bounds, 
in combination with Condition (f) in Proposition[2] will allow 
us to derive sufficient conditions for the existence of a Nash 
equilibrium in terms of structural properties of the network. 

F. The Stability of Nash Equilibria 

We present conditions for stability of a Nash equilibrium 
in terms of \ n (As). A Nash equilibrium x is asymptotically 
stable when the system of differential equations: 

/ii(x) = fi(x;5,A g ) 



■ Xl, 



x n = h n (x) = /„(x; S, Ag) - x n , 

is locally asymptotically stable around x. One can prove the 
following necessary and sufficient condition for an equilib- 
rium x to be asymptotically stable [5]: 

Lemma 5: An equilibrium profile x is asymptotically sta- 
ble if and only if 5 < — l/\ n (As) and #Xw=i dijXj > 1, 
for all inactive agents i £ Af\S. 

From the above lemma and Remark |4j we conclude that 
if 5 < — l/X n (Ag), there is a unique Nash equilibrium 
and it is asymptotically stable. These results show the close 
connection between the smallest and the largest eigenvalues 
of the adjacency matrix of interactions and the outcome of 



games with linear best response functions in a network. On 
the other hand, the results in this section are applicable if we 
are able to compute the eigenvalues of Ag. For many large- 
scale complex networks, the structure of the network can 
be very intricate [9]- [11] — in many cases not even known 
exactly — and an explicit eigenvalue decomposition can be 
very challenging to compute, if not impossible. In many 
practical settings, instead of having access to the complete 
network topology, we have access to local neighborhoods of 
the network structure around a set of nodes. In this context, 
it is important to understand the impact of local structural 
information on the eigenvalues of the adjacency matrix. In 
the rest of the paper, we propose a novel methodology to 
compute optimal bounds on relevant eigenvalues of Ag from 
local information regarding the network structure. Our results 
allow us to study the role of local structural information in 
the outcome of games with linear best response functions. 

III. Spectral Analysis of the Interaction Matrix 

We study the relationship between a network's local struc- 
tural properties and the smallest eigenvalue of its adjacency 
matrix. Algebraic graph theory provides us with tools to 
relate the eigenvalues of a network with its structural prop- 
erties. Particularly useful is the following well-known result 
relating the k-th spectral moment of Ag with the number of 
closed walks of length k in Q [7]: 

Lemma 6: Let Q be a simple graph. The fc-th spectral 
moment of the adjacency matrix of Q can be written as 



m k (Ac, 



1 



i 



(fc) 



(4) 



where is the set of all closed walks of length k in Q. 

From Q, we can easily compute the first three spectral 
moments of Ag in terms of the number of nodes, edges and 
triangles as follows [7]: 

Corollary 7: Let Q be a simple graph with adjacency 
matrix Ag. Denote by n, e and A the number of nodes, 
edges and triangles in Q, respectively. Then, 



mi{Ag) = 0, 
m 2 (Ag) = 2e/n, 
m 3 (A g ) = QA/n. 



(5) 



Remark 8: Notice that the coefficients 2 (resp. 6) in the 
above expressions corresponds to the number of closed walks 
of length 2 (resp. 3) enabled by the presence of an edge (resp. 
triangle). Similar expressions can be derived for higher-order 
spectral moments, although a more elaborated combinatorial 
analysis in required. 

In our case, we are interested in the following expressions, 
derived in [8], for the first five spectral moments of Q: 

Lemma 9: Let Q be a simple graph. Denote by e, A, Q 
and II the total number of edges, triangles, quadrangles and 

3 We denote by \Z\ the cardinality of a set Z. 



pentagons in Q, respectively. Define W% = XaLi d\ an ^ 
Cat = Yli diUt where t.- L is the number of triangles touching 
node i. Then, 

m 4 (Ag) = 1 [8Q + 2W 2 - e] , (6) 

n 

m 5 (A g ) = - [10n + WC dt - 30 A] . (7) 
n 

Observe how, as we increase the order of the moments, 
more complicated structural features appear in the expres- 
sions. In the following example, we illustrate how to use our 
expressions to compute the spectral moments of an online 
social network from empirical structural data. 

Example 10: In this example, we study a subgraph of 
Facebook obtained by exploring the online social network 
around a particular node, as follows. From a particular 
starting node, we crawl the network (using a breadth-first 
search) until we discover the set of all nodes that are within 
a radius 2 from the starting node. Using this set of nodes, 
and their interconnections (friendships), we construct a social 
subgraph T that has 2,404 nodes and 22, 786 edges. Using 
this real dataset, we compute the degrees di, the number 
of triangles t i7 quadrangles and pentagons pi touching 
each node i e V (F). In Fig. 1, we plot the distributions of 
degrees and triangles, as well as a scatter plot of U versus di 
(where each point has coordinates (di,ti), in log-log scale, 
for all i 6 V (J 7 ))- We can aggregate those quantities that 
are relevant to compute the spectral moments to obtain the 
following numerical values: 

e/n = J2 di/2n = 9.478, A/n = £ U/3n = 28.15, 
Q/n = £ ft/4" = 825.3, TJ/n = J2P*/5 n = 31, 794, 
W 2 /n = J2 d 2 Jn = 1, 318, C dt /n = £ d.tjn = 8, 520. 

Hence, using Corollary [7] and Lemma [9] we obtain the 
following values for the spectral moments: mi (Ajr) = 0, 
m 2 (Ajr) = 18.95, m 3 {Ajr) = 168.9, m 4 (A r ) = 9,230, 
and m 5 (Ajr) = 402,310. 

In this section, we have derived expressions to compute 
the first five spectral moment of Ag from network structural 
properties. In the next section, we use semidefinite program- 
ming to extract information regarding eigenvalues of interest 
from a sequence of spectral moments. 

IV. Optimal Spectral Bounds from Spectral 
Moments 

Here, we introduce an approach to derive an upper 
bound on the smallest eigenvalues of Ag from its sequence 
of spectral moments^] Since we have expressions for the 
spectral moments in terms of local structural properties, 
our bounds relate the eigenvalues of a network with these 
properties. There is a large literature studying the relationship 
between structural and spectral properties of graphs (see 
[12], [13], and references therein, for an extensive list of 

4 As a by-product of our analysis, we also derive lower bounds on the 
spectral radius of Ag, although these bounds are not essential in our 
analysis. 



spectral results). For many real-world networks, there is a 
particular set of structural properties that play a key role in 
the network's functionality. For example, it is well-known 
that social networks contain a large number of triangles 
(and other cycles). Hence, it would be useful to have 
spectral bounds where these structural features are jointly 
represented. In this section, we derive new upper bounds on 
the smallest eigenvalue of the adjacency matrix in terms of 
the structural properties involved in d3l, (|6]l and (|7]). Our 
results can be easily extended to derive lower bounds on the 
spectral radius of the adjacency matrix, although this bound 
is not of relevance in our analysis of games with linear best 
responses. 

Now, we derive bounds on the smallest eigenvalue of the 
adjacency matrix in terms of relevant structural properties 
by adapting the optimization framework proposed in [14]. 
We first need to introduce a probabilistic interpretation of a 
network eigenvalue spectrum and its spectral moments. For 
a simple graph Q, we define its spectral density as: 

1 ™ 

MaO = ~yX*-A<), (8) 

i=l 

where 8 (•) is the Dirac delta function and {A^}" =1 is the 
set of (real) eigenvalues of the (symmetric) adjacency matrix 
Ag. Consider a random variable X with probability density 
fig. The moments of X ~ fig are equal to the spectral 
moments of Ag, i.e., 

E M6 (X k ) - f x k fig Or) dx 

JR 

i ™ r 

= l x k S(x - Xi)dx 

n i=l 

1 ™ 

= ~J2*>i =m k (Ag), 
n *■ — ' 

i=l 

for all k > 0. 

In [14], Lasserre proposed a technique to compute the 
smallest interval [a, b] containing the supporj^] of a positive 
Borel measure fi from its complete sequence of moments 
(m r ) r>0 . In our spectral problem, the positive Borel measure 
under consideration is the spectral density fig (x), defined 
in p). Hence, in the context of our problem, the sequence 
of moments {m r ) r>0 is equal to (m r (Ag)) r>Q , and the 
smallest interval [a, b] containing the support of fig (x) is 
equal to [A„,Ai], by the definition in ([8|. 

Lasserre also proposed in [14] a numerical scheme to 
compute tight bounds on the values of a and b when 
a truncated sequence of moments (m r ) 0<r<k is known. 
This numerical scheme involves a series of semidefinite 
programs (SDP) in one variable. As we show below, at 
step s of this series of SDP's, we are given a sequence 
of moments (mi, TO2 S +i) and solve two SDP's whose 
solution provides an inner approximation [a s , j3 s ] Q [a, b]. In 
our case, since we have expressions for the first five spectral 

5 Recall that the support of a finite Borel measure fi on R, denoted by 
supp (fi), is the smallest closed set B such that fi (R\B) = 0. 
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Fig. 1. In the left and center figures, we plot the distributions of degrees and triangles of the social network under study (in log-log scale). In the right 
figure, we include a scatter plot where each point has coordinates (<i;,ti), in log-log scale, for all the nodes in the social graph. 



moments, (mi (Ag) , 7775 (Ag)), we can solve the first two 
steps of this series of SDP's to find inner approximations 
[a s , f3 s ] C [A n , Ax]. In other words, the solution to the SDP's 
provide us with the bounds a s > A n and f3 s < Ax. 

In order to formulate the series of SDP's proposed in [14], 
we need to define the so-called localizing matrix of our 
problem [15]. Given a sequence of moments, m( 2s+1 ) = 
(mi, m 2s+ i), our localizing matrix is a Hankel matrix 
defined as: 

H. (c) 4 R 2s+1 - c R 2s , (9) 

where R 2s and i?2s+x are the Hankel matrices of moments 
defined as 



1 

m 1 



TOx 

m 2 



in 



R 



2.S+1 — 



m s +i 
mi m 2 
m 2 m 3 



l s+l 



m 2s 



(10) 



m s+ i 
m s+ 2 



m s+ i m s+2 ■ ■ ■ m 2s+ i 

Hence, for a given sequence of moments, the entries of 
H s (c) depend affinely on the variable c. We can compute a s 
and f3 s as the solution to the following semidefinite programs 
[14]: 

Proposition 11: Let m' 2s+1 ' = (mi, m 2s+1 ) be the 
truncated sequence of moments of a positive Borel measure 
jit. Then, 

a < a s = max {a : H s (a) 0} , (11) 

a 

b > 0, 4 mm {/? : -F s 08) 0} , (12) 

for [a, 6] being the smallest interval containing supp(fx). 

Remark 12: Observe that a s and /3 S are the solutions to 
two SDP's in one variable, since the constraint H s (a) )p 
(resp. — H s ((3) )?= 0) indicates that the matrix H s (a) (resp. 
— H s (/?)) is positive semidefinite and this matrix has affine 
entries with respect to a (resp. (3). Hence, they can be 
efficiently computed using standard optimization software 



(for example, CVX [16]). As we increase s in Proposition 11 
more moments are involved in the SDP's, and the resulting 
bounds become tighter, i.e., a s +x < a s and /3 s +x > Ps- 
In the context of our spectral analysis, the Borel measure 



in Proposition 11 corresponds to the spectral density of 
a graph Q, and the smallest interval [a, b] corresponds to 
[A n (Ag) , p (Ag)]. Thus, Proposition [TT|provides an efficient 
numerical scheme to compute the bounds a s > X n (Ag) 
and /3 8 < p(Ag). When we are given a sequence of five 



spectral moments, we can solve the SDP's in (111 and ( 12 1 
analytically for s = 2. In this case, the localizing matrix is: 



fla(c) 



TOx — c 777.2 

777 2 — C777 ! 777 3 

7773 — C7772 m.4 



C777x 7773 — C7772 
C777 2 777-4 — C m 3 
Cm?, 777,5 — C7774 



(13) 



As we proved in Section III the spectral moments in the 
localizing matrix depend on the number of nodes, edges, 
cycles of length 3 to 5, the sum-of-squares of degrees W 2 , 
and the degree-triangle correlation Cdt- 

Furthermore, for s = 2, the optimal values a 2 and (3 2 can 
be analytically computed, as follows. First, note that H 2 (c) > 
(resp. —H 2 (c) >p 0) if and only if all the eigenvalues of H 2 
are nonnegative (resp. nonpositive). For a given sequence of 
five moments, the characteristic polynomial of H 2 (c) can be 
written as 

4> 2 (A) 4 det (XI - H 2 (c)) = A 3 + ft (c) \ 2 +p 2 (c) X+p 3 (c) , 

where pj (c) is a polynomial of degree j in the variable 
c (with coefficients depending on the moments). Thus, by 
Descartes' rule, all the eigenvalues of H 2 (c) are nonpositive 
if and only if pj (c) > 0, for j = 1,2, and 3. Similarly, 
all the eigenvalues are nonnegative if and only if p 2 (c) > 
and pi (c) ,P3 (c) < 0. In fact, one can prove that the optimal 



value of a 2 and j3 2 in ( 12 1 can be computed as the smallest 
and the largest roots of det H 2 (c) = 0, which yields a 
third degree polynomial in the variable c [14]. There are 
closed-form expressions for the roots of this polynomial (for 
example, Cardano's formula [17]), although the resulting 
expressions for the roots are rather complicated. 

In this subsection, we have presented a convex optimiza- 
tion framework to compute optimal bounds on the maximum 
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Fig. 2. Scatter plot of the spectral radius, p(d), versus the lower 
bound /?2 (Gi) , where each point is associated with one of the 100 social 
subgraphs considered in our experiments. 
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Fig. 3. Scatter plot of — A m i n (Gi) versus the bound —02 (Gi) , where 
each point is associated with one of the 100 social subgraphs considered in 
our experiments. 



and the minimum eigenvalues of a graph Q from a truncated 
sequence of its spectral moments. Since we have expressions 
for spectral moments in terms of local structural properties, 
these bounds relate the eigenvalues of a graph with its 
structural properties. 

V. Numerical Simulations 

As we illustrated in Section [II] there is a close connection 
between the largest and the smallest eigenvalues of a network 
and the outcome of a game with linear best response func- 
tions. In this section, we use our bounds on the support of the 
eigenvalue spectrum to study the role of structural properties 
in the existence and the stability of a Nash equilibrium. For 
this purpose, we analyze real data from a regional network 
of Facebook that spans 63, 731 users (nodes) connected by 
817, 090 friendships (edges) [18]. In order to corroborate our 
results in different network topologies, we extract multiple 
medium-size social subgraphs from the Facebook graph by 
running a Breath-First Search (BFS) around different starting 
nodes. Each BFS induces a social subgraph spanning all 
nodes 2 hops away from a starting node, as well as the edges 
connecting them. We use this approach to generate a set 
G = {Gi}i<ioo °f 100 different social subgraphs centered 
around 100 randomly chosen nodes H 

From Corollary [7] and Lemma |9j we can compute the 
first five spectral moments of a graph Gi from the fol- 
lowing structural properties: number of nodes (ni), edges 
(e,), triangles (A^), quadrangles (Qi), pentagons (II,), as 
well as the sum-of-squares of the degrees (Wi), and the 
degree-triangle correlation (Cj). For convenience, we define 
S(Gi) = {n^ei, Ai,Qi,Tli, Wi,C{\ as a set of relevant 
structural properties of Gi. In our numerical experiment, we 
first measure the set of relevant properties Si for each social 
subgraph Gi € G, and then compute the first five spectral 
moments of its adjacency matrix. From these moments, 

6 Although this procedure is common in studying large social network , 
it introduces biases that must be considered carefully [19]. 



we then compute the bounds a 2 (Gi) and /3 2 (Gi) using 



Proposition 11 As we mentioned before, these bounds can 
be computed as the maximum and minimum roots of a third 
order polynomial, for which closed form expressions are 
known. 

We illustrate the quality of our bounds in the following 
figures. Fig. 2 is a scatter plot where each circle has 
coordinates (p (Gi) , /3 2 (Gi)), for all Gi £ G. Observe 
how the spectral radii p (Gi) of these social subgraphs are 
remarkably close to the theoretical lower bound /3 2 (Gi). 
Therefore, we can use /3 2 (Gi) as an estimate of p(Gi) for 
social subgraphs. In Fig. 3 we include a scatter plot where 
each circle has coordinates (— \ n (Gi) , — a 2 (Gi)), for all 
Gi E G. Although a 2 (Gj) is a looser bound than /3 2 (Gi), 
we observe how there is a strong correlation between the 
value of A n (Gi) and a 2 (Gi). 

In these numerical experiments, we have first showed 
that «2 (Gi) and /? 2 (Gi) bound the smallest and the largest 
eigenvalues of the adjacency matrix, and that these bounds 
are tight, specially /3 2 (Gi). Since these bounds can be 
written as explicit functions of the structural properties in 
S (Gi), we can estimate the impact of structural perturbations 
on the spectral radius and the smallest eigenvalue by studying 
9/3 2 /5pi and <9a 2 / dpi for pi £ S(Gi). (Details of this 
perturbation analysis are left for future work due to space 
limitations.) 



VI. Conclusions 

In this paper, we have studied games with linear best 
response functions in a networked context. We have focused 
on analyzing the role of the network structure on the game 
outcome. In particular, the existence and the stability of a 
unique Nash equilibrium in this class of games are closely 
related with the smallest eigenvalue of the adjacency matrix 
of the network. We take this spectral result as the foundation 
to our work, and use algebraic graph theory and convex 



optimization to study how local structural properties of 
the network affect this eigenvalue. In particular, we have 
derived expressions for the first five spectral moments of 
the adjacency matrix in terms of local structural properties. 
These structural properties are: the number of nodes and 
edges, the number of cycles of length up to 5, the sum-of- 
squares of the degrees, and the degree-triangle correlation. 
From this sequence of five spectral moments, we propose 
a novel methodology to compute optimal bounds on the 
smallest and the largest eigenvalues of the adjacency matrix 
by solving two semidefinite programs. In our case, we are 
able to find analytical solutions to these optimal bounds by 
computing the roots of a cubic polynomial, for which closed- 
form expressions are available. Finally, we have verified the 
quality of our bounds by running numerical simulations in 
a set of 100 online social subgraphs. For future work, we 
shall use the results herein presented to study the effect 
of structural perturbations in the relevant eigenvalues of the 
adjacency matrix, and in properties of the Nash equilibrium. 
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